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I. INTRODUCTION 



Charmonium production at high energies has provided a useful laboratory for testing the 
high-energy limit of quantum chromodynamics (QCD) as well as the interplay of perturbative 
and non-perturbative phenomena in QCD. The factorization formalism of non-relativistic 

n 

QCD (NRQCD) |l| is a theoretical framework for the description of heavy- quarkonium 
production and decay. The factorization hypothesis of NRQCD assumes the separation of the 
effects of long and short distances in heavy-quarkonium production. NRQCD is organized 
as a perturbative expansion in two small parameters, the strong-coupling constant a s and 
the relative velocity v of the heavy quarks. 

The phenomenology of strong interactions at high energies exhibits a dominant role of 
gluon interactions in quarkonium production. In the conventional parton model |2j, the 
initial-state gluon dynamics is controlled by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equation [^. In this approach, it is assumed that S > fi 2 ^> Aq CD , 
where yS is the invariant collision energy, /i is the typical energy scale of the hard in- 
teraction, and Aqcd is the asymptotic scale parameter. In this way, the DGLAP evolution 
equation takes into account only one big logarithm, namely lnf^/Agco)- in fact, the collinear 
approximation is used, and the transverse momenta of the incoming gluons are neglected. 

In the high-energy limit, the contribution from the partonic subprocesses involving t- 
channel gluon exchanges to the total cross section can become dominant. The summation of 
the large logarithms lnf^v^S/ in the evolution equation can then be more important than 
the one of the ln(/i/AqcD) terms. In this case, the non-collinear gluon dynamics is described 
by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation "\. In the region under 
consideration, the transverse momenta (fc-r) of the incoming gluons and their off-shell proper- 
ties can no longer be neglected, and we deal with reggeized t-channel gluons. The theoretical 
framework for this kind of high-energy phenomenology is the so-called /^-factorization ap- 
proach j^, E|, which can be based on effective quantum field theory implemented with the 
non-abelian gauge-invariant action, as was suggested a few years ago Q. 

This paper is organized as follows. In Sec. |H1 the fc^-factorization approach is briefly 
reviewed and compared with the collinear parton model. The NRQCD formalism applied 
to heavy-quarkonium production is briefly recapitulated in Sec. IIHI In Sec. IIV1 we present 
in analytic form the squared amplitudes for S- and P-wave quarkonium production via the 
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fusion of reggeized gluons at leading order (LO) in ot s and v. In Sec. El we perform fits to the 
transverse-momentum (px) distributions of inclusive charmonium production measured at 
the Fermilab Tevatron to obtain numerical values for the non-perturbative matrix elements 
(NMEs) of the NRQCD factorization formalism. In Sees. EU and EEQ 

we compare our 

theoretical predictions with recent experimental data of charmonium production in 77, 7^, 
and deep-inelastic ep scattering at the DESY HERA and CERN LEP2 colliders. Section lVIIII 
contains our conclusions. 



II. THE FACTORIZATION APPROACH 

In the phenomenology of strong interactions at high energies, it is necessary to describe 
the QCD evolution of the gluon distribution functions of the colliding particles starting 
from some scale fio, which controls the non-perturbative regime, to the typical scale n 
of the hard-scattering processes, which is typically of the order of the transverse mass 



Mt = y M 2 + \p T \ 2 of the produced particle (or hadron jet) with (invariant) mass M 
and transverse two-momentum p^. In the region of very high energies, the typical ratio 
x = becomes very small, 1 <C 1. This leads to large logarithmic contributions of 

the type [a s ln(l/x)] n , which need to be resummed. This is conveniently done by adopting 
the high-energy factorization scheme, which also known as the /^-factorization approach, 
in which the incoming t-channel gluons have a finite transverse two-momentum and 
are off mass shell. This implies the notion of an unintegrated gluon distribution function 
$(x, |kr| 2 ,/i 2 ). The resummation is then implemented by the BFKL evolution equation 

Effective Feynman rules for processes involving incoming off-shell gluons were provided 
in Ref . |f| . The special trick is to choose the polarization four- vector of the incoming gluon 
as 

e\k T ) = *L y (1) 

where fc T = (0, k^, 0) is the transverse four- momentum of the gluon. In the case of gluon- 
gluon fusion, the four-momenta of the incoming gluons can be written as 

ki — X1P1 -\- k^rp^ 

k*2 = ^2-^2* "I" ^2Ti (2) 

where Pf = (VS/2)(1, 0, 0, 1) and P 2 M = (>/S/2)(l, 0, 0, — 1) are the four-momenta of the 



colliding protons in the center-of-mass frame. In the following, we shall also use the short- 
hand notation pr = |pr| etc. for the absolute of the transverse two- momentum. 

In Ref . p| , the incoming off-shell gluons are considered as Reggeons (or reggeized gluons) , 
which are interacting with quarks and on-shell Yang-Mills gluons in a specific way. Recently, 
in Ref. the Feynman rules for the effective field theory based on the non-abelian gauge- 
invariant action [7j were derived for the vertices RRg, Rgg, RRgg, Rggg, and RRggg, where 
R is an off-shell reggeized gluon and g is an on-shell Yang-Mills gluon. The interaction of a 
reggeized gluon with a quark is mediated via the transition vertex Rg. For the relevant LO 
a,„p.itudes, which are calculated below, both approaches flQ give the same answers. As 
was shown in Ref. [10j, the effective vertex RRg 8] can be obtained using the prescription 

for the off-shell gluon polarization four- vector of Eq. 

In the /^-factorization approach, which is based on the high-energy limit of QCD, the 
hadronic cross section of quarkonium (7i) production through the process 

p + p^h + x (3) 

and the partonic cross section for the reggeized-gluon fusion subprocess 

R + R^H + X (4) 

are related as 

da^(p + p^H + X,S)= J^-J d\k 1T \ 2 J Ikrrl 2 , /i 2 ) 

x /— /d|k 2r | 2 /^$(a; 2 ,|k 2r | 2 , / i 2 )dcT( J R + J R^H + X,k 1T ,k 2r ,S), (5) 

where s = XiX 2 S" — (kiy + k 2 ^) 2 , xi j2 are the fractions of the proton momenta passed on 
to the reggeized gluons, and y?i j2 are the angles enclosed between ki )2 T and the transverse 
momentum of 7i, which we take to point along the x axis. 

In our numerical calculations, we use the unintegrated gluon distribution functions by 
31umlein (JB) 11 1, by Jung and Salam (JS) Q], and by Kimber, Martin, and Ryskin (KMR) 



. A direct comparison between different unintegrated gluon distributions as functions 
of x, |k^| 2 , and fi 2 may be found in Ref. Note, that the JB version is based on the 

BFKL evolution equation jj| . On the contrary, the JS and KMR versions were obtained using 
the more complicated Catani-Ciafaloni-Fiorani-Marchesini (CCFM) evolution equation [15], 
which takes into account both large logarithms of the types ln(l/a;) and ln(yu/AQ C D)- 



For jj, ^> Aqcd and not too small x = fx/yS, the collinear approximation of the conven- 
tional parton model is recovered. In the collinear parton model, the hadronic cross section 
da{p + p — » TC + X, S) and the relevant partonic cross section dcr(g + g — > TC + X, s) are 
related as 

da PM (p + p^H + X,S) = Jd Xl G(x l7 /2 2 ) J dx 2 G(x 2 , ^ 2 )da(g + g -> ft + X, §), (6) 

where s = X\X 2 S and G(x,u 2 ) is the collinear gluon distribution function of the proton, 
which satisfies the DGLAP |3j evolution equation. The collinear and the unintegrated gluon 
distribution functions are formally related as 

xG{x^ 2 ) = rf|k T | 2 $(x, |k T | 2 ,/i 2 ), (7) 
Jo 

so that the normalizations of Eqs. (0) and © agree. 



III. NRQCD FORMALISM 

In the framework of the NRQCD factorization approach 1], the cross section of heavy- 
quarkonium production via a partonic subprocess a + b — > TC + X may be presented as a 
sum of terms in which the effects of long and short distances are factorized as 

da(a + b^H + X) =Y,da{a + b ^ QQ[n]+ X){O n [n]), (8) 

n 

where n denotes the set of color, spin, orbital and total angular momentum quantum numbers 
of the QQ pair and the four-momentum of the latter is assumed to be equal to the one of the 
physical quarkonium state Ti. The cross section da(a + b — ► QQ[n] + X) can be calculated 
in perturbative QCD as an expansion in a s using the non-relativistic approximation for 
the relative motion of the heavy quarks in the QQ pair. The non-perturbative transition 
of the QQ pair into 7i is described by the NMEs (O n [n\), which can be extracted from 
experimental data. 

To LO in v, we need to include the cc Fock states n = 3 S ( [\ 3 Sf\ 1 sf\ 3 P < § ,) if 
H = J/iprf and n = 3 P i j\ 3 Sf ) if H = XcJ, where J = 0,1,2. Their NMEs satisfy 
the multiplicity relations 

( iK«S)[3p(8)]) = (2J+ 1)(C^" 5 )[ 3 P (8) ]), 

[ 3 Pj 1} ]> = (2J+ 1)(C^°[ 3 P {1) ]), 

(0 X ^[ 3 S{ 8) }} = (2J+ 1){0^°[ 3 S{ 8) }}, (9) 
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which follow to LO in v from heavy-quark spin symmetry. For example, in the case of J /if) 
production, the wave function of the physical orthocharmonium state can be presented as a 
superposition of the Fock states: 

+ 0(v 2 )\cc[ 1 Sj 8) ]g) + ---, (10) 

where we use usual spectroscopic notation for the angular-momentum quantum numbers of 
the QQ pair and the index in parentheses (1, 8) denotes the color state, either color singlet 
or color octet. The color-singlet model (CSM) only takes into account the first term 
in Eq. (Ill))) , which is of order v°. In this case, the NME (O j ^[ 3 S[ 1] }) is directly related to 
the J /if) wave function at the origin \I/(0), which can be calculated in the framework of the 
quark potential model [17| . as 

(O-Wpsf 5 ]) = 2N C (2J+ 1)|^(0)| 2 , (11) 

where N c = 3 and J = 1. Similarly, the color-singlet P-wave NME reads 

(CF-fpW]) =2iV c (2J+l)|$ , (0)| 2 , (12) 



where ^/'(O) is the derivative of the Xcj wave function at the origin. 

In the general case, the partonic cross section of quarkonium production receives from 
the QQ Fock state n = 2S + 1 £ J y 8 ' the contribution 



of qua i 

flfi 



da(a + b - QQl^L^] - H) = da(a + b - QQ^L™]) ^ [ - ^ J) , (13) 

where A^ co i = 2N C for the color-singlet state, iV co i = — 1 for the color-octet state, and 
Np i = 2 J + 1. The partonic cross section of QQ production is defined as 

da(a + b - QQ[ 2S+1 L ( j fi) }) = j\A(a + b^ QQ[ 2S+1 L { j ,8) ])\ 2 d®, (14) 

where / is the flux factor of the incoming particles, which is taken as in the collinear parton 
model p (for example, I = 2X1X2S for process ©), A(a + b -> QQ[ 2S+1 L { }' 8) }) is the 
production amplitude, the bar indicates average (summation) over initial-state (final-state) 
spins and colors, and <i$ is the phase space volume of the outgoing particles. This convention 
implies that the cross section in the /^-factorization approach is normalized approximately 
to the cross section for on-shell gluons when ki^ = k 2 T = 0. 



j ]) can be obtained from the one for 



The production amplitude A(a + b — > QQ[ 2S+1 l[j 
an unspecified QQ state, ^4.(a + 6 — >■ QQ), by the application of appropriate projectors. The 
projectors on the spin- zero and spin-one states read l^ : 

1 /jp \ ( P 



n 



(/ - /// | --, \'- + q + m 



1 



Vt 



P ~ \ c 

q — m 7 



- + g + m 



(15) 



respectively, where p = 7 M p M , p M is the four-momentum of the pair, is the four- 
momentum of the relative motion, m = M/2 is the mass of the quark Q, and M is the 
mass of the quarkonium state Ti. In our numerical calculations, we use m c = 1.55 GeV. The 
projection operators for the color-singlet and color-octet states read: 

5i ; 



C-i 
C 8 -- 



(16) 



respectively, where T a with a = 1, . . . , N 2 — 1 are the generators of the color gauge group 
SU(3). To obtain the projection on a state with orbital-angular-momentum quantum number 
L, we need to take L times the derivative with respect to q and then put q = 0. For the 
processes discussed here, we have 



A(a + b 
A(a + b^QQ[ 3 S[ l ' 8) )) 

A(a + b^QQ[ 3 P^ 8) }) 



Tr 

Tr 

d 
dqp 



C h8 U A(a + b^QQ) 



g=0 



C lj8 U1A(a + b^QQ)e a (p) 



g=0 



■Tr 



C 1 , 8 U^A(a + b^QQ)e ocP (p) 



g=0 



(17) 



where e a (p) is the polarization four- vector of a spin-one particle with four-momentum p M 
and mass M = p 2 and £ a /3(p) is its counterpart for a spin-two particle. For the 3 Si state, 
the polarization sum reads 

PaPa' 



J2 £ a(p) £ *a>(p) = Vaa'ip) = "V + 



M 2 ' 



(18) 



For the 3 Pj states with J = 0, 1, 2, we have 



-(0) 



£4>)4?(p) = ^[v aa >{p)v^p) + v api {p)v al p{p)}-{v a pv al ^p). 



-(2)* 



(19) 



The subprocesses relevant for our analysis read: R + R — > QQ, R + R — > Q<5 + 
# + 7 -> <3<2, R + ^^QQ + g, R+e^e+ QQ, and R + e ^ e + QQ + g. 

IV. CHARMONIUM PRODUCTION BY REGGEIZED GLUONS 

In this section, we obtain the squared amplitudes for inclusive charmonium production 
via the fusion of two reggeized gluons or a reggeized gluon and a real or virtual photon in 
the framework of NRQCD. We work at LO in a s and v and consider the following partonic 
subprocesses: 

R + R - H[ 3 P?, 3 sf ] , l S { *\ 3 pf], (20) 
R + R - H[ 3 S { ^] +g, (21) 

J2 + 7 -> ft[ 3 S?VsfM 8) ], (22) 

i? + 7 - Wp^J+ff, (23) 

i? + e - e + Wpfif, 1 ^,^], (24) 

i? + e -> e + W[ 3 sf } ] + p. (25) 

Notice that, in the collinear parton model, subprocesses (|2()|). (|2*2*|l . and (|2*4*|) only contribute 
for w 0. Therefore, to LO in the collinear parton model, we need to take into account 
the corresponding subprocesses with an additional hard gluon in the final state, for example 
g+g — > H[ 3 S[ 8 ^]+g. The amplitudes of these color-octet subprocesses, after replacing g — » R 
in the initial state, are of next-to-leading order (NLO) in the /c^-factorization approach and 
suffer from infrared divergences, in contrast to the subprocesses (12 lj) and (|2*3*|) in the color- 
singlet channel. The analysis of NLO contributions to inclusive charmonium production by 
reggeized gluon-gluon fusion in the /c^-factorization approach is beyond the scope of this 
paper and needs a separate investigation. 



The phenomenological procedure, adopted in Ref . 20] , to regularize infrared divergences 
due to propagators getting on-shell with the help of some cut parameter, which is unknown 
a priori, is likely to be problematic. The analysis of NLO corrections in the /cy-factorization 



apprqac 
Ref 



F. [n| 



1 is currently an open issue, which has been consistently solved only in part, e.g. in 



where NLO corrections to the subprocess R + R — *> q were studied. 

n 

According to the prescription of Ref. p, the amplitude of R + R — * c + c(+g) is related 



S 



to the one of g + g — > c + c(+g) by 

A(R + R -> c + c(+<?)) = e^(ki)£ u (k 2 )A lJ , u (g + g^c + c(+g)), (26) 

where and e^{k 2 ) are defined according to Eq. ((TJ). Analogous relations hold for 

P + 7 — > c + c(+g) and R + e — > e + c + c(+g). The amplitudes of the relevant QCD 
subprocesses g + g — > c + c(+g) , g + 7 — > c + c(+(?), and g + e — > e + c + c(+g) are evaluated 
using the conventional Feynman rules of QCD. 

We now present and discuss our results for the squared amplitudes of subprocesses (I2UJ) 
and (J2U, contributing to hadropro duct ion. In the case of the 2 — > 1 subprocesses (J2Uj) . we 
obtain 



|.4(P + R -> W[ 3 P 



12 



2 2 



M 5 



|„4(P + P -> H[ 3 P 2 l 



(1)1 12 



|.4(P + P^ft[ 3 S{ 



(8)112 



|^(P+P-W[3P«]| 2 = fnW a { ° H ^ V ^fa.fr.y), 

32 2 ^(O^pP^j) r 3p] 
45 s J (^i^2,^), 

— Pi (t!,t 2 ,^), 



1 q( 8 )"||2 



|^(P + P^H[ 3 P (8) ]| 2 = 5vr 2 a 



.2 2 



|^(P + R -f ft[ 3 P^ 



(8)112 



|^(P + P -> ft[ 3 P 2 l 



(8)112 



2 2 
-7r a c 



(27) 



where 

P [35ll (*l,*2,¥>) 
P^l,^) 

P [3pol (ii,*2^) 

P [3Pll (^1^2,^) 



(M 2 + |p T | 2 ) [(ti + t 2 f + M 2 (*! + t 2 - 2 v ^t 2 "cos 



M 2 



(M 2 + tx + t 2 ) 2 
2 



(M 2 + *! + t 2 ) 



(M 2 + 



|Pt 



2\2 



Pt\ 1 sin y? 

(3M 2 + t! + t 2 ) cos <f + 2v^i2" 



2 M 2 (M 2 
9 

2 M 2 (M 2 + |p T | 2 ) 2 [(ti + t 2 ) 2 sin 2 <£> + M 2 (tj + 1 2 - 2 v ^t 2 "cos 
9 



(m 2 + ti + t 2 y 



1 



M 2 



3 (M 2 + tx + t 2 ) z 



(M 2 + |p T | 2 ) 2 {3M 4 + 3M 2 (tx + t 2 ) + 4txt 2 
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+ (*i + t 2 ) 2 cos 2 ip + [3M 2 + 2(ti + t 2 )] cos . (28) 

Here pt = kiT + k2T, £1,2 = |ki,2T| 2 , and 90 = y>i — ip 2 is the angle enclosed between kiT and 
k 2 T, so that 



|Pr| = *i + *2 + 2v^2"cos^. (29) 

It is interesting to consider the contribution of the diagram involving a three-gluon vertex 
separately. It is equal to 



S_ , 2<^[ 3 <Sf]> 



|^a(i2 + J2 - WpsT])! 2 = A, 2 ^^ 1 J/ (M 2 cos 2 v + | Pr | 2 ). (30) 



For |pt| 2 3> M 2 , one has 



\As(R + J2 -> HfSrW « ^ 2^3 Ipt|2 ' (31) 

which makes up the bulk of the contribution and can be interpreted as being due to the 
fragmentation production of the Ji meson. In fact, the right-hand side of Eq. (}3T| can be 
written in the factorized form 



\MR + H[ 3 S[ 8) W « \A(R + R^ g)\ 2 P(g - Hfs{ 8) }), (32) 

where 



\A(R+R^g)\ 2 = -na.\pr\ 2 (33) 

n 

refers to real-gluon production by reggeized-gluon fusion [10] and 

P(9^nM S) ]) = ^ a { ° H ^ )]) (34) 

is the probability for the fragmentation of a gluon to a 7i meson, which may be gleaned 
from the result for the corresponding fragmentation function at the starting scale //o j^j, 

D g^H[^]( z ^o) = ™ s — 5(1 - z). (35) 

The counterparts of Eq. ()27|) in the collinear parton model of QCD emerge through the 
operation 

\A(g + g^Hl^L^]\* = lira F ^ F ^\A(R + R - H^L™]\*. (36) 

<i,*2-»OJo 2,71 JO 2,71 
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In this way, we recover the well-known results 



\A(g + g^H[ 3 P, 



(i 



\A(g + g^H[ 3 Pi 



(i 



\A(g + g^H[ 3 P 



(i 



\A(g + g^H[^ 8 



\A(g + g^H[ 3 P, 



(8 



(8 



(8 



° 2 2 



0. 
32 
45 
0. 

5 



2 2 

7T a c 



2 2 

7r a: 



12 



tr 2 2 
07T a c 



0. 



(^[ 3 Pq (1) ]) 
M 3 

M 3 

M 

(^[ 3 Pq (8) ]) 
M 3 

[3p(8)l 



4 2 2 (Q"[ 3 P 2 W ]> 
3 71 °- M 3 



In the case of the 2 — > 2 subprocess (f2*Tj) . we find 



+ i2 -> H[ 3 S{ L> } + g\ 2 = n 3 a 
-320M 4 



, S <0 W N 1) ]> 



M 3 



81(M 2 - s) 2 (M 2 + ti - t) 2 (M 2 + * 2 



: (4t 1 t 2 M i (t 1 + t 2 + My cos 4 ((^i - y> 2 ) 



- 2 V / M 2 ~(ti + t 2 + M 2 ) cos 3 (^! - ^ 2 ){-t 2 2 M 2 - 3M 6 + tji + 3M 4 i + t 2 M 2 u + 3M A u 
+ i 2 tu - M 2 iu + t\{t 2 - M 2 + u) + hitl + i(M 2 + u)+ t 2 (5M 2 + i + u)) 

- 2 > /fc 2 "|p r |(f? - (t 2 - M 2 ){M 2 - t) + h{t 2 + 2M 2 + £)) cos^ 2 ) 

- (-M 2 + £ + uf(t 2 2 M 2 + t 2 M 4 + t\(t 2 + M 2 ) + t 2 t 2 + M 2 i 2 + M 2 u 2 

+ ^(t 2 + M 4 + u 2 + 2t 2 (4M 2 + £ + u)) - 2y/T 2 \v T \{tl + (t 2 + M 2 ){M 2 - t) 

+ h(t 2 + 8M 2 + t + 2u)) cos ^ 2 + (t 2 + (t 2 + M 2 ) 2 + 2*i(t 2 + 5M 2 ))|p T | 2 cos 2 <p 2 ) 

+ cos 2 ((^i - <^2)(* 2 M 4 - 2t 2 M 6 + 5t 2 M 8 - 2t 3 M 2 t + 4t 2 M 4 t - 10t 2 M 6 t + t 3 i 2 

- 2t\M 2 i 2 + 5t 2 M 4 t 2 - t 3 M 2 w - 7t 2 M 6 u + t\tu - t\M 2 iu + llt 2 M 4 £u - 3MHu 
+ t 2 2 Pu - At 2 M 2 i 2 u + 3M 4 P« + 2t 2 M 4 u 2 + t\iu 2 + 3M 4 M 2 - 2M 2 i 2 u 2 

+ t 3 (-4t 2 2 + (M 2 - u)(M 2 -t-u)+ t 2 (-5M 2 + t + u)) + tl(-4t 3 2 - 2M 6 + 4M 4 « 
+ iu(i + u)- 2t\ (12M 2 + i + u) - M 2 u(i + 2u) + t 2 (-25M 4 + t 2 + 7*u + 2m 2 
+ M 2 (8t + 15m))) + ti(5M 8 - 4M 2 M 2 + t 3 (-5M 2 + i + u) - M 6 (7t + 10m) 
+ M 4 (2t 2 + lltu + 5m 2 ) + t 2 (-25M 4 + 2t 2 + Itu + u 2 + M 2 (15t + 8u)) 
+ t 2 (-HM 6 + 15M 4 (t + u) + 6tu(t + u) - 2M 2 (4t 2 + 3tu + 4u 2 ))) 
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+ 2Vh\p T \(tt + t?(3t 2 + HM 2 + § + *) + * 2 (3* 2 + 7M 4 + si + * 2 (33M 2 - 2s + 3*) 
+ M 2 (lls + 4*)) + tx(tl - 3M 6 - 4M 2 * 2 + 7M 4 (s + *) + * 2 (5M 2 + s + 3*) 
+ * 2 (27M 4 + 2M 2 (s - 10t) + 6s*)) - (M 2 - *)(* 3 - t 2 M 4 + * 2 s 

+ M 2 (-4*(s + t) + M 2 (3s + At)))) cosy9 2 - 4* 2 |p T | 2 (* 3 + 2t\{t 2 + 2M 2 ) + 2(M 3 - Mtf 
+ + 2* 2 M 2 + 5M 4 - 4M 2 *)) cos 2 <p 2 ) - \p T \ 2 cos 2 ^i(2t% + At\t\ + 2t x t\ 
+ 4* 2 * 2 M 2 + 8*i* 2 M 2 + * 2 M 4 + 12*!* 2 M 4 + * 2 M 4 + 6*iM 6 + 10* 2 M 6 + M 8 
+ 4M 6 |p T | 2 - 2* 2 M 2 * - At x t 2 M 2 i - 2t 2 2 M 2 i - At x M 4 t - 20* 2 M 4 * - 2M 6 i 

- 8M 4 |p T | 2 * + t 2 t 2 + 2*1*2? + t 2 2 t 2 + 2*iM 2 * 2 + 10* 2 M 2 * 2 + M 4 * 2 + 4M 2 |p T | 2 * 2 

- 2* 2 M 2 w - 12*i* 2 M 2 w - 2* 2 M 2 m - 12*iM 4 w - 20* 2 M 4 m - 2M 6 u - 8M 4 |p T | 2 w 
+ 2t\iu + At x t 2 iu + 2t\iu + At x M 2 iu + 20* 2 M 2 *w + 2M 4 iu + 8M 2 \p T \ 2 iu + t\u 2 
+ 2*i* 2 m 2 + * 2 w 2 + 6hM 2 u 2 + m 2 M 2 u 2 + M 4 w 2 + 4M 2 |p T | 2 w 2 

- 8 v / tTM 2 |p T |(t 2 + M 2 — m)(M 2 - i - m)(cos</?i + cos(y?i - 2y? 2 )) 

+ 4 v / *I* 2 ~cos(y9i - y? 2 ) (M 2 - * - m) ((*! + * 2 ) 2 + M 2 (5M 2 + 2ti + 6* 2 - 4m)) 

+ 2*i cos(2(v?i - y? 2 ))(* 2 ((*i + t 2 ) 2 + M 2 (2tx + 4* 2 + 5M 2 - 4m)) + 2M 2 (m 2 - uf) 

- 16M 2 |p T |v^cosv? 2 (M 2 -£-M) 2 +4M 2 |p T | 2 cos(2v? 2 )(M 2 (M 2 -2t-2M) + (* + w) 2 ) 

- 2 v / tT(M 2 - * - u) cos(v?i - y? 2 )(-2|p T |(2* 2 * 2 — M 6 + 2* 2 (M 2 - *) + M 4 i+M 4 u - 2M 2 tu 
+ * 2 (8M 4 - 7M 2 * - *(* + u)) + *i(2* 2 + M 2 (3M 2 - * - 3m) + * 2 (11M 2 + * + 3m))) cos y? 2 

+ ^ft 2 (2t\M 2 + 7* 2 M 4 — M 6 + * 2 |p T | 2 + 2* 2 M 2 |p T | 2 + 5M 4 |p T | 2 

+ * 2 (2* 2 + 2M 2 + |p T | 2 ) - 4* 2 M 2 * - 4M 2 |p T | 2 * + * 2 * 2 + 3M 2 * 2 - t 2 M 2 u - t 2 iu 

- M 2 tu - t 2 u + 3MV - iu 2 + *i(2* 2 + 7M 4 + M 2 (6|p T | 2 - i - Au) + u(-i + u) 
+ * 2 (13M 2 + 2|p T | 2 + 3(* + «))) + |p T | 2 (t 2 + t\ + 2t 2 M 2 + 5M 4 + 2*i(* 2 + 3M 2 ) 

- 4M 2 t) cos(2^ 2 ))) + 2|pr| cos^i(2*iV^(*i + t 2 + M 2 )(t\ + M 4 — M 2 u 

+ *i(* 2 — M 2 +u) + t 2 {2M 2 + u)) cos 3 (v?i - v? 2 ) + v / ^iCos 2 (v?i - <^ 2 )(* 2 + ll^M 2 + 7^M 4 

- 3t 2 M 6 + t\s + ll* 2 M 2 s + 7* 2 M 4 s - 3M 6 s + t\u + 4* 2 M 2 m + 7* 2 M 4 m - AM 6 u 
+ t 2 2 su + 7M 4 sm - 4* 2 M 2 m 2 + 8M 4 i 2 - 4M 2 sm 2 - AM 2 u 3 + t\{t 2 - M 2 + u) 

+ * 2 (3* 2 + s(-M 2 +u)+ * 2 (5M 2 + s + 3m)) + *i(3* 3 + M 6 — M 4 m 

+ * 2 (33M 2 - 2s + 3m) + * 2 (27M 4 + 2M 2 (s - 10m) + 6su)) + 2^/t 2 ~\p T \(-5M 6 + 5M 4 * 

- 4* 2 (M 4 - M 2 t) + t\{M 2 -i-u) + t\(M 2 -t-u) + 5M 4 m - 2t 1 (2M 2 (M 2 - u) 
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+ t 2 (M 2 + i + u))) cos <f 2 ) + Vhi-M 2 + i + u) 2 (tit 2 + t\ + tiM 2 + 8t 2 M 2 + M 4 

+ 4M 2 |p T | 2 + 2t 2 i - hu + t 2 u - M 2 u - 2 v / £ 2 ~|pt|(7M 2 + i + u) cos^ 2 + 4M 2 |p T | 2 cos(2y? 2 )) 

- (M 2 — i — u) cos(^i - y? 2 )(-(|p T |(3M 6 -3M 4 £- 12t 2 (M 4 - M 2 £) +t 2 (M 2 - t - u) 
+ t 2 (M 2 - t - «) - 3M 4 u + 4M 2 M - 2ti(6M 2 (M 2 - u) + t 2 (HM 2 + 3(i + u)))) cos<^ 2 ) 

- 2y/t 2 {2t\{t 2 + M 2 -u)- M 2 (M 4 + 2t{\p T \ 2 + u)- M 2 (2|p T | 2 + 1 + u) 
+ t 2 (-3M 2 + 3t + u)) + ti(2t 2 + 8M 4 + M 2 (2|p T | 2 - 7u) - u{t + u) 

+ t 2 (llM 2 + 3t + u)) + 2M 2 | Pr | 2 (t 1 + M 2 — *) cos(2</? 2 ))))), (38) 

where s = (k± + k 2 ) 2 , i = (ki —p) 2 , and u = (k 2 —p) 2 are the standard Mandelstam variables. 
With the aid of Eq. ()36j) . we recover from Eq. the well-known collinear-parton- model 



result 



23] 



\ A ( g + g ^ n[ 3 S U}+g\* = ttV 



(O h [ 3 S[ 1] ]) 320M 4 



M 3 81(M 2 -t) 2 (M 2 -w) 2 (t + m) 2 
x (M 4 £ 2 - 2M 2 t 3 + t 4 + M 4 £m - 3M 2 t 2 u + 2t 3 u + M 4 u 2 

- 3M 2 iu 2 + 3t 2 u 2 - 2M 2 u 3 + 2tu 3 + u A ). (39) 



We now turn to subprocesses (|22|) and (j23|) . with one real photon in the initial state. For 
the 2 — > 1 subprocesses (J22J), which are pure color-octet processes, we find 

0, 

2 a (OTfif]> 

57r aa s en — , 

M 

32 2 2 (^[ 3 P (8) ]) (3M 2 + ^) 2 
3"^ M 3 (M 2 + t 1 ) 2 ' 
64 2 2 (gT^]) t!(2M 2 + t!) 
3 n aas6 Q M 3 (M 2 + t0 2 ' 

64 2 2 (C W [ 3 P 2 {8) ]) 6M 4 + 6M 2 tx + t 2 , N 

— n 2 aa s e 2 a - 1 2 J/ — ^ — i, 40 

15 M 3 (M 2 + ti) 2 ' v ; 

where eg is electric charge of the heavy quark Q. Ap plic ation of Eq. (|3l)j) to Eq. (J4T?|) yields 



I-40R + 7- 


^H[ 3 ^ 8) ]| 2 = 


\A(R + -f- 




\A(R + j- 


+ H[ 3 Pi 8) ]\ 2 = 


\A(R + j- 


+ H[ 3 P} S) ]\ 2 = 


\A(R + j- 


->H[ 3 P 2 {8) }\ 2 = 



the well-known results of the collinear parton model 



2i, 



|^ + 7 ^ W[ 3^)]|2 = o 



\A{g + j ^H[ 3 Pi S) }\ 2 = 96n 2 aa s e 2 Q 



M 3 
13 



128 2 2 ((^[ 3 P 2 (8) ]) 



For the 2 — > 2 subprocess (J2*3|) . which is a color-singlet process, we find 

(O^S^]) 2048M 2 



M 3 27(M 2 - s) 2 (M 2 - w) 2 (t! + M 2 — t) 2 
x (t 4 M 2 + M 2 (s 2 + su + u 2 - M 2 (s + u)f + t 3 (M 2 (5s + 3u) - 7M 4 - su) 

+ t 2 (sw(w - s) + M 4 (3w - lis) + M 2 (7s 2 + 2su - 3m 2 )) + M(sw 2 + M A (u - 6s) 

+ M 2 (4s 2 + sm - m 2 )) - 2v^|p T |(t?M 2 + t 2 (-7M 4 - su + M 2 (3s + 4m)) 

+ ti(M 4 (-7s + 2u) - s 2 m + M 2 (2s 2 + su - 2m 2 )) - M 2 (2M 4 (s + u) - 2M 2 m(3s + 2u) 

+ u(3s 2 + 4su + 2u 2 ))) cos^ - 2M 2 |p r | 2 (t 3 + M 2 s 2 + t\(M 2 + 2s) 

+ ti(2M 2 s + s 2 -2t 2 ))cos 2 ^), (42) 



where k 2 now represents the photon four-momentum and (p 2 = 0. Equation (|42|) agrees with 



the corresponding result in Ref. 25[, but has a more compact form. By means of Eq. (|36|h 



Eq. (|42jl collapses to the well-known collinear-parton model result 

(O^pSf 5 ]) 2048M 4 



\ A{g + 1 ^ H[ 3 S m ]+g \2 = n 3 aa 2 e 2 Q 



M 3 27(M 2 -t) 2 (Af 2 -u) 2 (t + u) 2 
x (M 4 ? - 2M 2 t 3 + t 4 + MHu - 3M 2 i 2 u + 2t 3 u + M A u 2 

- 3M 2 iu 2 + 3t 2 u 2 - 2M 2 u 3 + 2tu 3 + u 4 ). (43) 

Finally, we turn to subprocesses ()24j) and (|2"5)l . through which electroproduction proceeds 
at LO. As for the 2^2 subprocesses (|24j) . which are all color-octet processes, we have 



\A{R + e-+e + H[ 3 S{ 8) }\ 2 = 0, 
lAiR + e^e + H^S^l 2 = 64vr 3 « 2 a s e 2 1> 



m 3 y 2 Q 2 (M 2 + Q 2 + t 1 

x ((2 + (y a - 2)y 2 )((M 2 + t,) 2 + Q 4 + 2Q 2 M 2 - 2Q 2 t iy2 ) + AQ% 



+4QyJt 1 (l - y 2 )(M 2 + Q 2 + t x ){y 2 - 2)y 2 cos ( Vl - <p 2 ) 

+ 2{y 2 - 1) (Q 4 + (M 2 + t x ) 2 + Q 2 (2M 2 + 2t x - t lV 2 2 )) cos (2(^ - <^ 2 ))W 



|.4(P + e->e + ft[ 3 P (8) ]| 2 = ^7r 3 a 2 a s e 



256_ 3 , - <(^[ 3 P (8) ]> 1 



3 Q M 5 y 2 Q 2 (M 2 + Q 2 + ^) 4 

x ( f 2 + ( !, 2 - 2)!, 2 ) ( 9.1 r , 2 \M\Q 2 + h) + 22M 4 Q 4 + 22M 4 t 2 + (Q 2 + tx) 2 (Q 4 + t 2 ] 
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+ 8M 2 (Q 2 + ti)(Q 4 + t 2 )) + 2M 4 g 2 tx(52 + y 2 ((43 - 9y 2 )y 2 - 64)) 
+ 2Q 2 t x (g 2 + t0 2 (10 - y 2 (U + (y 2 - 6)y 2 )) 
+ AM 2 Q 2 t 1 (Q 2 + ix)(16 - 3y 2 (8 + (y 2 - % 2 )) 

+ 4Qy / t 1 (l-y 2 )(M 2 + g 2 + ti)(g 4 (4 + (y 2 - 2)y 2 ) + 2g 2 (t x (4 + (y 2 - % 2 ) 

+ M 2 (8 + y 2 (-4 + 3y 2 ))) + (3M 2 + i 1 )(i 1 (4 + (y 2 - 2)y 2 ) 

+ M 2 (4 + y 2 (3y 2 - 2)))) cos - y? 2 ) - 2(3M 2 + Q 2 + - 1) 

x (3M 6 + 7M 4 (g 2 + h) + (g 2 + ti)(g 4 + 1\ + g 2 tx(2 + (y 2 - 4)y 2 )) 

+ M 2 (5g 4 + U\ + Q 2 *i(10 + y 2 {3y 2 - 4)))) cos (2(^ - v? 2 )))m 2 , 

^_ 512 ^ 2 . , 2 (^[ 3 A (8) ]) i 



x ((2 + (y 2 - 2)y 2 ) (Q s + h(M 2 + t 1 ) 2 (2M 2 + h)) 
+ 2g 6 (y 2 - 2)(M 2 (y 2 -2)-t 1 (2-y 2 + yl)) 

+ g 4 (4M 2 t 1 (y 2 - 3)(y 2 - 2) + M 4 (10 + (y 2 - W)y 2 ) - 2t 2 (y 2 (6 + y 2 (-5 + 2y 2 )) - 6)) 
+ 2g 2 (M 4 t 1 (10 + (-8 + y 2 )y 2 ) - 2M\y 2 - 1) - ^(y 2 - 2)(2 + (y 2 - l)y 2 ) 

- M 2 tl(y 2 (10 + y 2 (2y 2 - 5)) - 12)) + AQ^l - y 2 )(M 2 + Q 2 + t 1 ){M\y 2 - 2) 
+ (g 2 + t 1 ) 2 {y 2 - 2)y 2 - M 2 (g 2 (2 + y 2 ) + t 1 (2 + y 2 - 2y 2 2 ))) cos - </? 2 ) 

+ 2(y 2 - i)((g 2 + to 4 - g 2 t x (g 2 + to 2 ^ 2 + m\(q 2 + ^) 2 - 2g 2 t 11/2 ) + 2M 2 ((g 2 + 

- Q\(Q 2 + h)y 2 - Q 2 t\yD) cos (2( Vl - p 2 ))) M 2 , 

3p(8)l 



512 , ((^pP^]) 



+ e ^ e + HM*\\ 2 = ^«W 2 A^_^ ^(M^g^tO 4 
x ((2 + (j/ 2 - 2)y 2 ) (g 8 + (M 2 + ti) 2 (6M 4 + QM 2 h + t 2 )) + 2g 6 (M 2 (8 + (y 2 - 8)y 2 ) 

- h(y 2 - 3)(y 2 - 2) 2 ) + g 4 (M 4 (38 + y 2 (7y 2 - 38)) + 4M 2 ti(20 + y 2 (8y 2 - 25)) 

+ t 2 (44 - 2y 2 (30 + y 2 (2y 2 - 13)))) + 2Q 2 {-{t\{y 2 - 3)(y 2 - 2) 2 ) + 6M 6 (3 + (y 2 - 3)y 2 ) 

- M 2 t 2 (y 2 (50 + y 2 (6y 2 - 25)) - 40) - M 4 t 1 ( 2/2 (52 + y 2 (6y 2 - 25)) - 46)) 

+ 4g v /t 1 (i-y 2 )(M 2 + g 2 + ti)(g 4 (4 + ( y2 - 2) V2 ) 

+ t 2 (4 + (y 2 - 2)y 2 ) + 3M 4 (2 + y 2 {2y 2 - 3)) + M 2 tx(10 + y 2 (6y 2 - 11)) 
+ g 2 (M 2 (10 - lly 2 ) + 2*x(4 + (j/ 2 - % 2 ))) cos - y? 2 ) 

- 2(y 2 - l)(2M 2 ((g 2 + tO 3 " 5g%(g 2 + h)y 2 + 3Q 2 tlyl) 



15 



+ (Q 2 + ti) 2 (Q 4 + t\ + Q\(2 + (y 2 - 4)2/3)) 

+ M 4 (Q 4 + t\ + 2Q\(1 + 3(2/ 2 - l)y 2 ))) cos (2(^ - </? 2 )))m 2 . (44) 

As usual, Q 2 = — q 2 and y 2 = (q ■ P) / (k ■ P), where P M , A/\ and q^ = — are the 
four-momenta of the incoming proton, the incoming lepton, the outgoing lepton, and the 
virtual photon, respectively, tpi is the angle between ki^ and py, and (p 2 is the angle between 
qr and py. The corresponding formulas in the collinear parton model J27] are recovered as 
explained in Eq. ()36|) and read: 



\A(g + e - 


^e + H[ 3 s[ 8) }\ 2 


\A{g + e - 


+ e + W[ 1 SJ 8) ]| 2 


\A(g + e- 


->e + H[ 3 P (8) ]| 2 


\A(g + e- 


+ e + H[ 3 Pi 8) ]\ 2 



\A(g + e^e + HliP 2 {8) ]\ 2 = ^nWa s e 2 Q 



Our analytic result for the 2 
refrain from listing it here. 



0. 



64:7r 3 a 2 a s eQ 



(O^S®]) y 2 -2y 2 



M 



y\Q 2 



256 



(O h [ 3 P, 



3p(8)j 



M 3 



— 7T a a s e Q - 
(y 2 -2y 2 + 2)(Q 2 + 3M 2 ) 2 



y 2 Q 2 (Q 2 + M 2 ) 2 
3" 71 ° as6Q M 3 

2 o„. 1 oA/02 /i/„, i\N/\#2 



512 



X 



((y 2 -2y 2 + 2)Q 2 - A(y 2 - 1))M 2 



512 



y 2 2 (Q 2 + M 2 ) 2 
(C[ 3 P 2 (8) ]> 



M 3 



((y 2 - 2y 2 + 2)(Q 4 + 6M 4 ) - 12(y 2 - 1)M 2 Q 2 ) 

ylQ 2 {Q 2 + M2 ) 2 ' 

3 color-singlet subprocess (|25|) is rather lengthy, and we 



V. CHARMONIUM PRODUCTION AT THE TEVATRON 



During the last decade, the CDF Collaboration at the Tevatron 28LI29I ] collected data on 
charmonium production at energies v^S* = 1.8 TeV (run I) and \fS = 1.96 TeV (run II) in 
the central region of pseudorapidity |r/| < 0.6. The data cover a large interval in transverse 
momentum, namely 5 < < 20 GeV (run I) and < pr < 20 GeV (run II). The data 
sample of run I [28] includes pr distributions of J / ip mesons that were produced directly in 
the hard interaction, via radiative decays of Xcj mesons, via decays of ip' mesons, and via 
decays of b hadrons. That of run II 2^] includes px distributions of prompt J / ip mesons, so 
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far without separation into direct, Xcj-decay, and ip' -decay contributions, and of J/ip mesons 
from 6-hadron decays. 

As is well known, the cross section of charmonium production measured at the Tevatron 
is more than one order of magnitude larger than the prediction of the CSM evaluated 
within the collinear parton mode 



/cT-factorization approach (20 , 



31 



]30(. Switching from the collinear parton model to the 
somewhat ameliorates the situation, but still does not 
lead to agreement at all. On the other hand, a successful description of the data could be 
achieved with the NRQCD factorization formalism implemented in the collinear parton 



model, including the fusion and fragmentation mechanisms of charmonium hadroproduction 



133 



34 ]. 



Charmonium hadroproduction was studied some time ago using the NRQCD factorization 



formalism implemented in th e /c^ -factorization approach invoking both the fusion 
and fragmentation pictures 



te kr- 



20 



31 



that, in order to describe the 



It was found l2C 

experimental data from the CDF Collaboration |2j|, it is necessary to employ a set of NMEs 
that greatly differs from the one favored by the collinear parton model. In this paper, we 
confirm this conclusion only to some degree. 

On the other hand, the polarization of prompt J/ip mesons measured at the Tevatron 



35| also provides a sensitive probe of the NRQCD mechanism. This issue was careful! 



investigated both in the collinear parton model 



3q and in the /c^-factorization approach (37 1. 



None of these studies was able to prove or disprove the NRQCD factorization hypothesis. 

In contrast to previous analyses in the collinear parton model or the /cr-factorization 
approach, we perform a joint fit to the run-I and run-II CDF data 3,13 ^° obtain the color- 
octet NMEs for J/ip, ip', and Xcj mesons. We use three different versions of unintegrated 
gluon distribution function. Our calculations are based on exact analytical expressions for 
the relevant squared amplitudes, which were previously unknown in literature. Our fits 
include five experimental data sets, which come as pt distributions of J/ip mesons from 
direct production, prompt production, Xcj decays, and ip' decays in run I and from prompt 
production in run II. 

We now describe how to evaluate the differential hadronic cross section from Eq. (JSJ in 
combination with the squared matrix elements of the 2 — ^ 1 and 2 — > 2 subprocesses (I20J) 
and (|2ip. respectively. The rapidity and pseudorapidity of a charmonium state with four- 
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momentum p^ = (p°, Pt,P 3 ) are given by 

In 



y 



1 . p° + p 3 



2 p° - p 3 ' 



1 p \+p 3 
- m 
2 



Pi -P 3 ' 



(46) 



respectively. For the 2 — > 1 subprocess (J20J), we have 



da KT (p + p-+ ft + X) 


Ipt 


d|p T |cfy 


(|p T | 2 + M 2 ) 2 



IT 



#1 



x $ p (6, IkxTl 2 ,^)^^, |k 2T | 2 , / i 2 )|^l( J R + i? - W)P 



where 



6 



p° + p 3 



6 



(47) 



(48) 



and k 2 r = Pt — k^. in our numerical analysis, we choose the factorization scale to be 
/i = Mt- For the 2 — > 2 subprocess (j2Tf . we have 

rfa KT (p + p^H + X) 



d|p r |dy 



|Pr| 

(2vr) ; 



d\ki T \ 2 J dcpi J dx 2 J d\k 2 T\ 2 J d(f 2 



x %( Xl , |k 1T | 2 , /i 2 )%(a; 2 , |k 2T | 2 , ^ 



where 



1 



Xi 



(x 2 - f 2 ) (2x1X25) ; 



^ki T + k 2 T - Pt) 2 - M 2 - |p T | 2 + a^iS 1 



(49) 



(50) 



We now present and discuss our results. In Table HI we list out fit results for the relevant 
color-octet NMEs for three different choices of unintegrated gluon distribution function, 
namely JB Q, JS Q, and KMR . The color-singlet NMEs are not fitted, but deter- 
mined from the measured partial decay widths of ip(nS) — > l + + l~ and % c2 ^7 + 7. 
The numerical values are adopted from Ref. Q and read: (O^S^}) = 1.3 GeV 3 , 
(O^Fsf 5 ]) = 6.5 x KT 1 GeV 3 , and {0^[ 3 Pf ] }) = (2 J + 1) x 8.9 x 1(T 2 GeV 5 . They 
were obtained using the vacuum saturation approximation and heavy-quark spin symm etry 
in the NRQCD factorization formulas and inclu ding NLO QCD radiative corrections 
The relevant branching ratios are taken from Ref. 
B(i\)' -> J/ip + X) = 0.576, B( X co -> J/?P + 7) = 0.012, B( Xc i -> J/ip + 7) = 0.318, and 
B(xc2 —> J/tp + 7) = 0.203. They somewhat differ from the values used previously j^J. For 



0.0601, 



comparison, we list in Table U also the NMEs obtained in Ref. 34j for the collinear parton 
model with the LO parton distribution functions of the proton by Martin, Roberts, Stirling, 
and Thorne (MRST98LO) |4l|. 
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We first study the relative importance of the different intermediate states in direct J/ if) 



28|, with p T > 5 GeV, the 



and if)' production. In previous fits to CDF data from run I 
linear combinations 

M? = (O^S^}) + -^(O n [ 3 P^}) (51) 

m c 

for 7i = J /if), if)' were fixed because it was infeasible to separate the contributions propor- 
tional to (O^S^}) and <O w [ 3 P {8) ]). By contrast, the new run-II data which reach 
down to pt = 0, allow us to determine (O n [ 1 S^}} and (O h [ 3 Pq]) separately because 
the respective contributions exhibit different pt dependences for px < 5 GeV. This feature 
is nicely illustrated in Fig. ^ where the shapes of the relevant color-octet contributions to 
prompt J /if) production, proportional to (O n [ 3 S\ 8) ]), (O^S^]), and (C W [ 3 P (8) ]), are com- 
pared with that of the CDF data from run II 29] . Notice that the color-octet contributions 
differ in the peak position, by up to 1 GeV. Apparently, this suffices to disentangle the con- 
tributions previously combined by Eq. flEE). We find that (£> J/W/ [ 3 P (8) ]) and (O^'^S^]) 
are compatible with zero, independent of the choice of unintegrated gluon density — a strik- 
ing result. For the case of J /if) production from if)' decay, this implies that the 3 Si and 
3 S[ 8 ^ channels are sufficient to describe the measured px distribution (see Fig. 0J). 

In Figs.EHSl we compare the CDF data on J /if) mesons from direct production, if)' decays, 
and Xcj decays in run I -28] and from prompt production in run II j^J, respectively, with the 
theoretical results evaluated with the NMEs listed in Table U] From Fig. El we observe that 
the color-singlet contribution is significant, especially at low values of pt, and comparable 
to the one from the 1 S , q 8 ' ) channel. As is familiar from the collinear parton model, the 3 S[ 8 ^ 
contribution makes up the bulk of the cross section at large values of pt- Incidentally, the 
values of (O j ^[ 3 S[ 8 ' 1 }) obtained in the /cy-factorization framework are in average quite close 
to the one obtained in the collinear parton model, as may be seen from Table HI The situation 
is very similar for J /if) production from if)' decay, considered in Fig. El except that the 
and 3 P < f ) contributions are negligible. 

At this point, we wish to compare our results for direct J /if) hadropro duct ion in the 
fcT-factorization approach with the literature, specifically with Refs. j^, |^ , which consider 
the partonic subprocess (J20J) - By contrast, in Ref. [^jj, the NLO subprocess R + R — > 
J/ip^Si^] +g was studied, leaving aside the LO subprocess (pH) . In Ref. j^], the value 
(O j ^[ 3 S{ 8) }) = 7.0 x 1(T 3 GeV 3 was obtained using the Kwiecinski-Martin-Stasto (KMS) 



42j unintegrated gluon distribution function. This value is 2.6 times larger than the result 
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we found using the KMR version, which is very similar to the KMS one. We attribute 
this difference in (0 J ^[ 3 S[ 8 ^}} to the different scale choice, /x = kx, used by the authors 
of Ref. Adopting their value for (0 J ^[ 3 S (8) }), we can reproduce their result for the 

respective cross section contribution. On the other hand, the value (0 J ^[ 3 S[ 8S) }} = 15.0 x 
10~ 3 GeV 3 found in Ref. bol exceeds the one of Ref. \M by a factor of 2.1 and our KMR 

n 

value by a factor of 5.6. Furthermore, the cross section evaluated in Ref. |20J falls off with 
Pt considerably more slowly than in Ref. [&| and here, only by one order of magnitude as 
Pt runs from 2 to 20 GeV, while the unintegrated gluon density in the proton falls off with 
kr far more rapidly. 

The discussion of J / if) production from radiative Xcj decays, considered in Fig. HI is sim- 
pler because there is only one free parameter in the fit, namely (O Xc0 [ 3 S[ 8 ^]}. We confirm the 
conclusion of Ref. that, in the fey-factorization approach, the color-singlet contribution 
is sufficient to describe the data. In fact, the best fit is realized when (O Xe0 ^Si ]) is taken 
to be zero or very small. In case of the JB gluon density, the fitting procedure even favors 
a negative value of (O Xc0 [ 3 S{ 8) ]). 

In Fig. 03 the px distribution of prompt J / if) production in run II is broken down into the 
contributions from direct production, if)' decays, and Xcj decays. We observe that the latter 
is dominant for ~ 5 GeV, while prompt J / if> mesons are preferably produced directly at 
larger values of px- The contribution from if)' decays stays at the level of several percent 
for all values of px- While the JS [3] and KMR gluon densities allow for a faithful 
description of the measured px distribution 3] , the JB [□] one has a problem in the low-py 
range, at px ^ 5 GeV, where even the Xcj-decay contribution, which is entirely of color- 
singlet origin, exceeds the data. This problem can be traced to the speed of growth of the 
JB gluon density as kx — > 0. By contrast, the JS and KMR gluon densities are smaller and 
approximately kx independent at low values of kx- For this reason, we excluded the CDF 
prompt- J /if) data from run I j3] and run II [3] from our fit based on the JB gluon density. 

Considering the color-octet NMEs relevant for the J/ip, if)' and Xcj production mecha- 
nisms, we can formulate the following heuristic rule for favoured transitions from color-octet 
to color-singlet states: AL ~ and AS* ~ 0; i.e. these transitions are doubly chromoelectric 
and preserve the orbital angular momentum and the spin of the heavy-quark bound state. 



20 



VI. CHARMONIUM PRODUCTION AT HERA 



At HERA, the cross section of prompt J/if> production was measured in a wide range of 
the kinematic variables W 2 = (P + q) 2 , Q 2 = —q 2 , y 2 = (P-q)/(P-k), z = (P ■ p)/(P ■ q), p T 
and y, where P M , g M = k^ — k' 11 , andp M are the four-momenta of the incoming proton, 

incoming lepton, scattered lepton, virtual photon, and produced J /if) meson, respectively, 
both in photoproduction at small values of Q 2 , and deep-inelastic scattering (DIS) [44], 



at large values of Q 2 . At sufficiently large values of Q 2 , the virtual photon behaves like a 
point-like object, while, at low values of Q 2 , it can either act as a point-like object (direct 
photoproduction) or interact via its quark and gluon content (resolved photoproduction). 
Resolved photoproduction is only important at low values of z. 

In the region z ^ 1, diffract ive production, which is beyond the scope of this paper, takes 
place. In order to suppress the diffractive-production contribution, one usually applies the 
acceptance cut z < 0.9. This effectively eliminates the contributions from the 2 —>■ 1 partonic 
subprocesses (|2*2*j) and (f24j) . so that we are left with the 2^2 partonic subprocesses (|2*3*j) 
and (125J) . 

Let us first present the relevant formulas for the double differential cross sections of DIS, 
direct photoproduction, and resolved photoproduction. In the case of DIS, we have 
da KT (p + e -> e + H + X) 1 



d\p T \ 2 dz 8z(2%) 



dQ 2 J dy 2 J d\k 1T \ 2 J d<pi J dtp 2 



where 



x 1 



x $ (r Ik I 2 ^ lAR + e^e + H + gW 

X <P p (£i, k 1T , \l ) — — — , [bz) 

(2/2 - X2){2xiS) 2 



[k 1T + q 2T - Pt) 2 - M 2 - \p T \ 2 + y 2 X\S + {V2 ~ X2)Q' 



(j/2 - X2)S 



Xi = 9F , X2 = ■ (53) 

ZH/p Ztli e 

Here, E p and E e are the proton and lepton energies in the laboratory frame, and we have 
S = AE p E e and |q 2T | = J(l - y 2 )Q 2 . 



In the case of direct photoproduction, we have 
da KT (p + e -> e + H + X) 



d\p T \ 2 dz 2z(2tt) 



dy 2 J d\k 1T \ 2 J dipi 



t • „ ,2 ju / AA(R + j^H + g 2 

x $ p (xi, k 1T 2 ,/x 2 )f l/e (y 2 )— f Vo ' » 54 

2/2(2/2 - X2){ZX 1 S) 2 
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where 



{V2 - X2)S 



[k 1T - p T ) 2 — M 2 — |p T | 2 + 2/2X15 



(55) 



and f 1 / e {y 2 ) is the quasi- real photon flux. In the Weizacker- Williams approximation, the 
latter takes the form 



a 



Z7T 



1 



1 



)2 

' max / 



(56) 



) 2 

' max 



IP- Qmin \ Qmin Qi 

where QL- = m 2 yl/ (1 — y 2 ) and Q max is determined by the experimental set-up, e.g. Q\ 
1 GeV 2 y . 

In the case of resolved photoproduction, we take into account the 2 — > 1 and 2-^2 
partonic subprocesses (|2Ti|) and (|2*T|). respectively, where the first reggeized gluon comes 
from the proton and the second one from the photon. For subprocess (|20p. the relevant 
doubly differential cross section reads: 

da KT (p + e -> e + ft + X) 1 



cft/ 2 / d\k 1T \ 2 J 



d\p T \ 2 dz 2z(\p T \ 2 + M 2 ) 2 

x%{ Xl , \k 1T \ 2 ,fx 2 )f j/e (y 2 )^(x 2 , \k 2T \ 2 ,fx 2 )\A(R + R^H)\ 2 , 



where 



«i = Xi, x 2 



X2 

y2 



^2T 



Pt - k 



IT- 



(57) 



(58) 



For subprocess 1)21)1. the relevant doubly differential cross section is given by 
da KT (p + e-> e + H + X) 1 



d\p T \ 2 dz 



2z(l-z)(2 



it 



dy 2 J d\k 1T \ 2 J d(fi J dx 2 J d\k 2T \ 2 J d(p 2 



x$ p (xi, |k 1T | 2 ,//)/ 7/e (2/ 2 )$ 7 (x 2 , |k 2T | 2 ,/i 2 )' 



x 2 (2x!X 2 y 2 Sy 



where 



£1 



'k iT - p T ) 



M 



|Pt| 2 + x 2 y 2 XiS 



(59) 



(60) 



(^22/2 - Xz)S 

To evaluate the unintegrated gluon distribution function in the resolved photon, 
$ 7 (x 2 , \k 2T \ 2 , fi 2 ), we use a procedure suggested by Blumlein j^j], which is similar to the 
proton case [11]. As input for this, we use the collinear parton distribution functions of the 
resolved photon by Gliick, Reya, and Vogt (GRV 7 ) 

In Figs. our NRQCD predictions in the /cy-factorization approach, evaluated with 
the NMEs from Table HI are compared with the HERA data (43, 44]. Specifically, Figs. El and 
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[7|refer to the p\ and z distributions in photoproduction with E p = 820 GeV, E e = 27.5 GeV, 



60 GeV < W < 240 GeV, and Q 2 < 1 GeV 2 



43j |. while Figs.|K]and|n]refer to those in DIS with 



K, = 920 GeV, E e = 27.5 GeV, 50 GeV < W < 225 GeV, and 2 GeV 2 < Q 2 < 100 GeV 2 



44j. Acceptance cuts common to both photoproduction and DIS include Pt > 1 GeV and 



0.3 < z < 0.9. In this regime, the LO NRQCD predictions in the /^-factorization approach 
are mainly due to the color-singlet channels and are thus fairly independent of the color- 
octet NMEs presented in Table U] Therefore, our results agree well with previous calculations 
in the CSM j^, up to minor differences in the choice of the color-singlet NMEs and the 
c- quark mass. 



VII. CHARMONIUM PRODUCTION AT LEP2 



Some time ago, the DELPHI Collaboration presented data on the inclusive cross section 
of J /if) photoproduction in 77 collisions (e + + e~ — > e + + e~ + J /if) + X) at LEP2, taken 
as a function of the J/ip transverse momentum pt j3- The J/ip mesons were identified 
through their decays to fi + fi~ pairs, and events where the system X contains a prompt 
photon were suppressed by requiring that at least four charged tracks were reconstructed. 
The average e + e~ center-of-mass energy was y/S = 197 GeV, the scattered positrons and 
electrons were antitagged, with maximum angle (9 max = 32 mrad, and the maximum 77 
center-of-mass energy was chosen to be W = 35 GeV in order to reject the major part of 
the non-two-photon events. 

Under LEP2 experimental conditions, most J j ip mesons are produced promptly, while the 
cross section for J I ifr mesons from 6-hadron decays is estimated to be about 1% of the total 
J j ip cross section |49J and can be safely neglected. Because the average value of the photon 
virtuality Q 2 is small, the Weizsacker- Williams approximation can be used to evaluate the 
e + e _ cross section from the 77 cross section as 

da(e + + e- -> e+ + e~ + H + X) = J dy x J dy 2 / 7/e ( 2 / 1 )/ 7/e (y 2 )d ( r(7 + 7 - H + X). (61) 

The process e + + e~ — >• e + + e~ + J/ip + X receives contributions from direct, single- 
resolved, and double-resolved photoproduction. The relevant partonic subprocesses are: 7 + 

7 - ^n+s, 1+ r - h?e$MK i+R - n 3 s?}+g, r+r - n 3 sf\ ^f, 3 p ( h 

and R + R — > H^S^] + g. The squared amplitude of 7 + 7 — > H[ 3 S[ 8 ^} + g may be found 
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in Ref . 49] , the ones for the other partonic subprocesses were presented in Sec. HV 
The cross section of direct photoproduction is evaluated as 
da(e + + e~ -> e+ + e~ +H + X) 1 



d\p T \ 2 dy 



4tt 



dy2f 1 / e (yi)f J /e(y2 



ym \A(r/ + 1 -^n + g) \ 2 



y2 - 6 (ZymS) 2 

where 6 and 6 are defined in Eq. ()48p and 

_ y&S - M 2 

1 (y2 - 6)5 ' 

In the case of single-resolved photoproduction via the 2 — > 1 subprocesses, we have 

da KT (e + + e~ -> e+ + e" + H + X) 
d\p T \ 2 dy 



(62) 



(63) 



4?r / dy 1 f l/e {y l )f l/e {y 2 



x$ 7 (xi, |k 1T | 2 ,/i 2 )2/ 2 



|.4(i? + 7 -> ft)| 2 



(64) 



(2x 1 y 1 y 2 S) 2 

where X\ = £1/2/1, 2/2 = £2, and k 1T = p T . In the case of single-resolved photoproduction 
via the 2 — > 2 subprocess, we have 

da KT (e+ + e- ^ e+ + e~ + ft + X) 1 



d\p T \ 2 dy 



2(2vr) 



2 , dyi I dy 2 I d\k 1T \ 2 \ d^ f l/e (y 1 )f l/e (y 2 ) 



x $Jx 1} |ki T | ,/x 



2 2 , 2/2 |.4(i? + 7 ~> ?* + <?) |< 



2/2 - 6 



where 



1 



(2x 1 y 1 y 2 Sy 



(kiT - Pt) 2 - M 2 - | Pr | 2 + y 2 65 



2/1(2/2 -6)5 

In the case of double-resolved photoproduction via the 2 — ^ 1 subprocesses, we have 

da KT {e + + e~ -> e+ + e" + W + X) 
d\p T \ 2 dy 



(65) 



(66) 



2 / dyi / C?2/2 / d|ki T | / d^l /-y/e (3/l ) /-y/e (s/2 ) 



x ^Jk^.^^lM'.^ ^^gl' , (67) 

where Xi = 6/2/I) x 2 = 6/2/2, and k 2 T = Pt — k ir . In the case of double-resolved photo- 
production via the 2^2 subprocess, we have 

da KT (e + + e- ^ e + + e- +H + X) 1 



d\ PT \ 2 dy 2 (2n)' J J J ^ f 

x J dx 2 J d\k 2T \ 2 J dLp 2 f^ /e (y 1 )^ y (x 1 ,\k 1T \ 2 ,fj 2 )f^ f/e (y 2 )^^(x 2 ,\k 2T \ 2 ,ij, 2 ) 

2/2 \A(R + R-+ H + g)\ 2 



x 2 y 2 - 6 (2x 1 x 2 y 1 y 2 S) 2 



(68) 
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where 

xi = — — \(k 1T + k 2T - Pr) 2 - M 2 - \p T \ 2 + x 2 y2^iS\ . (69) 

In Fig. we confront the p\ distribution of e + + e~ — > e + + e~ + J/ip + X, where X 
is devoid of prompt photons, measured by DELPHI with our full theoretical prediction 
(line No. 4), which is broken down into the single- resolved color-octet contribution (line 
No. 1), the single- resolved color-singlet contribution (line No. 2), and the direct plus double- 
resolved contributions (line No. 3). We observe that the single-resolved contribution makes 
up the bulk of the cross section, while the direct and double-resolved contributions are greatly 
suppressed, and that, within the single- resolved contribution, the color-singlet channel is 
dominant. The experimental data overshoot the theoretical prediction by a moderate factor 
of 2-3. For the case of 77 collisions, we conclude that the color-singlet processes are dominant 
in the /c^-factorization approach, a situation familiar from photo- and electroproduction in 
ep collisions considered in Sec. IVII The situation is quite different for the collinear parton 



model, where color-octet^ processes dominate 
Recently, in Ref. 



251 ] . it was attempted to interpret the DELPHI data in the kr- 
factorization approach invoking only the CSM and neglecting the cascade decays of the 
ip' and Xcj mesons. Curve No. 2 in Fig. El approximately agrees with the corresponding 
predictions in Ref. 25] for m c = 1.55 GeV. In Ref. 2^]i a significantly lower value of m c is 
employed to reach agreement with the DELPHI data. 

VIII. CONCLUSION 

Working at LO in the /^-factorization approach to NRQCD, we analytically evaluated 
the squared amplitudes of prompt charmonium production by reggeized gluons in RR, R'y, 
and Re collisions. We extracted the relevant color-octet NMEs, (C w [ 3 ^f } ]), (O^Sf } ]), and 
(C W [ 3 P (8) ]) for H = J/ip, V', and Xcj through fits to p T distributions measured by the CDF 



Collaboration in pp collisions at the Tevatron with y/S = 1.8 TeV |2j| and 1.96 TeV [29j 



using three different versions of unintegrated gluon distribution function, namely JB [ll( , JS 
[l^ . and KMR Appealing to the assumed NRQCD factorization, we used the NMEs 



thus obtained to predict various cross section distributions of prompt J/ ip photoproduction 
and electroproduction in ep collisions and photoproduction in e + e~ collisions and compared 
them with ZEUS Q and HI Q data from HERA and DELPHI H data from LEP2, 
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respectively. In the case of photoproduction, we included both the direct and resolved 
contributions. As for the unintegrated parton distribution functions of the proton and the 
resolved photon, we assumed the gluon content to be dominant. 

Our fits to the Tevatron data turned out to be satisfactory, except for the one to the 
Xcj sample based on the JB gluon density in the proton, where the fit result significantly 
exceeded the measured cross section in the small-pr region. We found agreement with the 
HERA and LEP2 data within a factor of 2, which is the typical size of the theoretical 
uncertainty due to the lack of knowledge of the precise value of the c-quark mass and 
the NLO corrections. Specifically, we found that direct and resolved photoproduction in 
ep collisions under HERA kinematic conditions dominantly proceed through color-singlet 
processes, namely R(p) + 7 -> H[ 3 S[ 1] ] + g and R(p) + R(<y) -> H[ 3 S[ 1} ) + g, respectively. 
Similarly, photoproduction in e + e~ collisions under LEP2 kinematic conditions is mainly 
mediated via the color-singlet subprocess R(^y) + 7 — > W[ 3 5i ] + g, but the color-octet 
subprocess R(j) +7 — > H[ Sq ] also contributes appreciably. 

LO predictions in both the collinear parton model and the /^-factorization framework 
suffer from sizeable theoretical uncertainties, which are largely due to unphysical-scale de- 
pendences. Substantial improvement can only be achieved by performing full NLO anal- 
yses. While the stage for the NLO NRQCD treatment of 2 — > 2 processes has been set 
in the collinear parton model j^J, conceptual issues still remain to be clarified in the kx- 
factorization approach. Since, at NLO, incoming partons can gain a finite kr kick through 
the perturbative emission of partons, one expects that essential features produced by the 
/cT-factorization approach at LO will thus automatically show up at NLO in the collinear 
parton model. 
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FIG. 1: Contributions to the px distribution of prompt J ftp hadroproduction in pp scattering with 
\/~S = 1.96 TeV and \y\ < 0.6 from the relevant color-octet states compared with CDF data from 
Tevatron run II 29]. All distributions are normalized to unity at their peaks. 
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FIG. 2: Contributions to the pt distribution of direct J /if) hadroproduction in pp scattering with 
VS = 1.8 TeV and \y\ < 0.6 from the partonic subprocesses (1) R + R -► J/ip[ 3 s[ 8) ], (2) it! + 
R -> J/i/jV-sfpfpjp], (3)_R + i? -c J/^[ 3 5i 1} ] + 5, and (4) their sum compared with CDF 



data from Tevatron run I 28 1. The theoretical results are obtained with the (a) JB [ill], (b) JS 

n n 

[12j, or (c) KMR [lj] unintegrated gluon distribution functions. The decay branching fraction 
B(J/tp — > /i + + is included. 
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FIG. 3: Contributions to the px distribution of J /ip mesons from tp' decays in hadroproduction in pp 
scattering with =1.8 TeV and \y\ < 0.6 from the partonic subprocesses (1) R + R — > ip'[ 3 Si], 
(2) R+R -> V't 1 ^ 8 ^ 3 ^?] ( this contribution actually vanished), (3) R+R -► ?//[ 3 ,sf and ( 4 ) 
their sum compared with CDF data from Tevatron run I 28]. The theoretical results are obtained 

n n n 

1 If ] - (b) JS [12], or (c) KMR [13j] unintegrated gluon distribution functions. The 



with the (a) JB 



decay branching fraction B(J/tp — » n + + /z ) is included. 
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FIG. 4: Contributions to the pt distribution of J /ip mesons from \cj decays in hadroproduction 
in pp scattering with = 1.8 TeV and \y\ < 0.6 from the sum of the partonic subprocesses 
R+R — ► Xcj[ 3 Pj ] an d R+R — > XcjP^ 8 ^], the latter of which being quite unimportant, compared 
with CDF data from Tevatron run I |28J|. The theoretical results are obtained with the JB ll il l . 

□ n 

JS [121 ] . or KMR [13( unintegrated gluon distribution functions. The decay branching fraction 
B{J/ij) — > /i + + is included. 
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FIG. 5: Contributions to the px distribution of prompt J ftp hadroproduction in pp scattering with 
■</S = 1.96 TeV and \y\ < 0.6 from (1) direct production, (2) ip' decays, (3) Xcj decays, and (4) 
their sum compared with CDF data from Tevatron run II 29]. The theoretical results are obtained 
with the (a) JB 11], (b) JS 12], or (c) KMR unintegrated gluon distribution functions. The 
decay branching fraction B(J/ip — > fi + + n~) is included. 
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FIG. 6: Contribution to the distribution of prompt J/ip photoproduction in ep scattering with 

E p = 820 GeV, E e = 27.5 GeV, 60 GeV < W < 240 GeV, Q 2 < 1 GeV 2 , and 0.3 < z < 0.9 
rom the direct-photon subprocess R + 7 — ► 7i[ 3 S^] + g compared with ZEUS data from HERA 
3- The resolved-photon subprocesses R + R^ WP fif > , 3 p[ 1] , 3 sf V S^L 3 pj 8) ] are neglected. 

The theoretical results are obtained with the JB [lj, JS [12j, or KMR [13J unintegrated gluon 

distribution functions. 
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FIG. 7: Contributions to the z distribution of prompt J/tp photoproduction in ep scattering with 
E p = 820 GeV, E e = 27.5 GeV, 60 GeV < W < 240 GeV, Q 2 < 1 GeV 2 , and (a) p T > 1 GeV, 
(b) pt > 2 GeV, or (c) px > 3 GeV from (1) the direct-photon subprocess R + 7 — ► 'H[ 3 S'| 1 ^ + g, 
(2) the resolved-photon subprocesses R + R^ H[ 3 s{ 1] , 3 pj 1} , 3 S^, 1 S^ 8 \ 3 pf], and (3) their sum 
compared with ZEUS data from HERA [2] . The theoretical results are obtained with the JB 11 1 
unintegrated gluon distribution function. 
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FIG. 8: Contribution to the pj, distribution of prompt J/tp electroproduction in ep scattering with 
E p = 920 GeV, E e = 27. .5 GeV, 50 GeV < W < 225 GeV, 2 GeV 2 < Q 2 < 100 GeV 2 , and 
0.3 < z < 0.9 from the color-singlet subprocess R + e — ► e + H[ 3 s[ ] + q compared with HI data 
from HERA 44]. The theoretical results are obtained with the JB 
unintegrated gluon distribution functions. 



-f q compared with HI data 
Q, JS Q, or KMR Q 
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FIG. 9: Contribution to the z distribution of prompt J /ip electroproduction in ep scattering with 
E p = 920 GeV, E e = 27.5 GeV, 50 GeV < W < 225 GeV, 2 GeV 2 < Q 2 < 100 GeV 2 , and 
px > 1 GeV from the color-singlet subprocess R + e — > e + W[ 3 5^] + o compared with HI data 
from HERA Q- The theoretical results are obtained with the JB |ll|, JS Q, or KMR Q 
unintegrated gluon distribution functions. 
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FIG. 10: Contributions to thep^ distribution of prompt J/ip photoproduction in e + e annihilation 
with y/S = 197 GeV, Q 2 < 9.93 GeV 2 , W < 35 GeV, and \y\ < 2 from the partonic subprocesses 
(1) R + 7 - H^sf ^pf], ( 2 ) R + 7 - + 5, (3) 7 + 7 - WfSf ] + 5, A + i* - 

W[ 3 5S 8) , 1 5S 8) , 3 P2lj and R + R ~* n[ 3 S { i } } + 9, and (4) their sum compared with DELPHI 
data from LEP2 |48|. The theoretical results are obtained with the JB [llj| unintegrated gluon 
distribution function. 
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TABLE I: NMEs for J/ifi, tp', and Xcj mesons from fits in the collinear parton model (PM) 
using the MRST98LO par ton distribution functions of the proton | and in the fc^-factorization 
approach using the JbQ, JsQ, and KM R ll3H unintegrated gluon distribution functions. The 
CDF prompt data from run I [28( and run II [29J have been excluded from our fit based on the JB 
gluon density. 



NME 



PM [34] 



Fit JB 



Fit JS 



Fit KMR 

1.3 
2.7 x 1(T 3 
1.4 x 1(T 2 


1.4 x 10 -2 



{0 J ^[ 3 s{ 1] ])/GeV 3 
(0 J ^[ 3 s[ S) ])/GeV 3 
(O^S^/GeV 3 
(O J /^[ 3 P (8) ])/GeV 5 



M 3 J/ //GeV 3 



1.3 

4.4 x 1(T 3 



1.7 x 10" 2 



1.3 

1.5 x 10~ 3 

6.6 x 1(T 3 



6.6 x 1(T 3 



1.3 
6.1 x 1(T 3 
9.0 x 10~ 3 


9.0 x 10~ 3 



(O^psfty/GeV 3 
(^'[ 3 5{ 8) ])/GeV 3 
(<9^'[isJ 8 )])/GeV 3 
(OV>'[3p ( 8 )])/GeV 5 
MfjGeV 3 



6.5 x 10" 1 
4.2 x 10~ 3 



1.3 x 10" 



6.5 x 10" 1 
3.0 x 10~ 4 









6.5 x lO" 1 
1.5 x 10~ 3 









6.5 x lO" 1 
8.3 x 10~ 4 









(0x<*[ 3 P o (1) ])/GeV 5 
(0^[ 3 S[ 8) ])/Ge\ 3 



8.9 x 10~ 2 
2.3 x 10~ 3 



8.9 x 10" 2 




8.9 x 10~ 2 
2.2 x 10~ 4 



8.9 x 10~ 2 
4.7 x 10~ 5 



X 2 /d.o.f 



2.2 



4.1 



3.0 
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